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Abstract 

Multi-system interaction is an important and difficult problem in physics. Motivated by the ex- 
perimental result of an electronic circuit element "Fractor" , we introduce the concept of dynamic- 
order fractional dynamic system, in which the differential-order of a fractional dynamic system is 
determined by the output signal of another dynamic system. The new concept offers a comprehen- 
sive explanation of physical mechanism of multi-system interaction. The properties and potential 
applications of dynamic-order fractional dynamic systems are further explored with the analysis of 
anomalous relaxation and diffusion processes. 

PACS numbers: 66.10.C-, 05.10.Gg, 02.60.Cb 
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Fractional dynamic system has been focused by physicists and mathematicians over the 
last decades, and has received great success in the analysis of anomalous diffusion [l-4], 



viscoelastic rheology js, 6], control systems |7|, |8j, complex networ 
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human tissue and electrochemical corrosion process 



lfj, etc. 
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cs [9 1, wave dissipation in 
13|. 



In the past several years, various physical applications have given birth to the variable- 



order fractional dynamic system 14H16j. However, the multi-system interaction must be 



considered in the physical mechanism analysis of the variable-order fractional dynamic sys- 
tem and its applications. Especially, the behavior of a dynamic system may change with the 
evolution of other dynamic systems in multi-system physical processes. How to characterize 
the interaction effect between these dynamic systems? Establishing a system of equations 
which includes intricate interaction terms, will cause great difficulties in modeling and com- 
putation. Even worse, since they may miss capturing the critical physical mechanism of 
the considered problems, the established model will produce incorrect results which greatly 
deviate from experimental results or field measurement data. Meanwhile, researchers have 
confirmed that the differential orders in some fractional dynamic systems are non-constant 



and are often functions of other variables or system outputs |16l . Il7j . For instance, Glockle 
and Nonnenmacher have found that the differential order of proteins relaxation is a function 
of temperature [181 ] . Therefore, in order to exploit the physical mechanism of a variable- 
order fractional dynamic system, another dynamic system usually should be considered. For 
simplicity, we name this type of variable-order fractional dynamic system as dynamic-order 
fractional dynamic system. 

The purpose of this letter is to make an innovative study of dynamic-order fractional dy- 
namic systems. Since previous studies have indicated that the differential-order is a critical 
factor in the variable-order fractional dynamic system, it is necessary to analyze the physi- 
cal mechanism about how environmental factors or system variables influence the system's 
differential order. From the analysis of fractional dynamic system, especially the variable- 



order dynamic system, it has been found the fractiona 
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differential order of one dynamic 
1 

18l |. Hence, in many fractional sys- 



system usually stems from another dynamic system 
terns, the differential order should be called dynamic-order. The comprehensive study of 
dynamic-order system will give us new understanding about multi-system interaction. 

Firstly, we introduce the motivation of dynamic-order fractional system via an experiment 
of Fractor. This experiment will give us a basic understanding of dynamic-order fractional 
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dynamic systems. 

A Fractor is a two lead passive electronic circuit element similar to a resistor or capacitor, 
exhibiting a non-integer order power- law impedance versus frequency 19|. The Fractor is 
proven effective as a feedback element in control systems for real-world applications such as 
temperature control, robotics, etc. 20(. The prototype Fractor is made by hand and not 
much larger than typical through-hole capacitors and the typical unit is 3.5 cm on a side 
and about 1.0 cm thick, as shown in Fig. [I] (Left). The impedance behavior of the Fractor 
can be accurately modeled by the following formula which is achieved by Laplace transform 
of the fractional order operator Zp rac (<jj) = K/(jcor) x , K is the impedance magnitude at 
calibration frequency (u c = 1/t); A G (0, 1) is the fractional exponent and u is the frequency 



19|, 



21] . From the previous study of the Fractor, it has been confirmed that the order of 



the Fractor may change with the temperature and its internal material structure. Several 
experimental results have implied that, temperature can influence the derivative order or 
integral order which will determine the capacity of the equipment 18|, ll9] . In our experiment, 
we intend to give a clear relationship between the temperature and the fractional order of the 
Fractor. In this experiment, temperature is controlled by the equipment of Quanser Heatflow 
Experiment (HFE) [19|. The evolution process of temperature is a typical heat transfer 
process. In this heat transfer system, the instantaneous temperature value of the Fractor is 
the output signal we intend to obtain. Based on the previous experimental observation, we 
rewrite the above formula as 

K 

Z FracH ~ J—pjxpy I 1 ) 

where A(T) e (0, 1) is the fractional exponent, T is the temperature. 

In the experiment, we firstly placed the Fractor in the HFE unit, then the temperature 
of the Fractor was measured by three temperature sensors. Next, the order of the Fractor 
was calculated by an HP DSA and fitted via the expression ([I]). Finally, The relationship 
between the temperature and the order of the Fractor was obtained as shown in Fig. [1] 
(Right). The experimental details were stated in 20] . It is observed that the order of 
the Fractor is an approximately linear function of the temperature. Since the order of the 
Fractor changes with the temperature which is the output signal of the heat transfer system 
in the Fractor, this system is a typical dynamic-order system. 

Next, we discuss the dynamic-order fractional dynamic system from the viewpoint of 
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FIG. 1: Left: Photo of a sample hand-made prototype Fractor. The detailed preparation and 
description of the Fractor can be found in [19|]. Right: The fractional order evolution data of 
the Fractor. In this temperature range, the linear fitting function is A(T) = p\T + P2, where 
pi = 0.001127 ± 0.000116, p 2 = 0.8835 ± 0.0048. 

variable-order fractional dynamic system. A representative definition of variable-order frac- 
tional derivative in Caputo sense can be stated as follows 22|, [23] 

1 ? /'(r)dr 



c D«*>f(t) 



< a{Z) < 1, 



(2) 



r(l-a(Z))7 (t-T)°W 
where Z denotes certain system variable or other independent variable, which govern the 
differential-order of interested dynamic-order system. at(Z) is a function of the indepen- 
dent variable Z and T() is the Gamma function. In some variable-order fractional dynamic 
systems, the variable-order is a function of certain variable, such as temperature or con- 
centration, while the order-influencing variable usually can be regarded as output signal of 
another dynamic system. 

In engineering practice, we usually encounter the phenomena of the multi-system interac- 
tion, such as temperature field, stress field and electromagnetic field coupled together with 
the concerned dynamic system under investigation. From the traditional approach, when 
considering the coupling effect of different systems or some multi-scale physical processes, 
we are accustomed to employing a system of differential equations and they have received 



great success in the past decades 



24]. For example, when modeling the tracer transport in a 



fractured granite, Reimus et al. |25| established a convection-diffusion equation in fracture 
and a diffusion equation in matrix, in which the intricate interaction is reflected by an ad- 
ditional term in the convection-diffusion equation in fracture scale. However, in the recent 
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FIG. 2: An illustration of a dynamic-order fractional dynamic system. In this schematic, u\ is the 
input of system 1 and at\(t) is the differential order, y\(t) is the output signal. 02(2/1 (t)) denotes 
the differential-order of system 2, u 2 and y% (t) represent the input and output signals of system 2. 

decades, experimental and theoretical studies have implied that the dynamic-order system 
approach can be employed to reflect the interaction effect in some real-world multi-system 



problems 16|, llSj. The link between these systems may be established via differential order 
of each dynamic system, which can be visually illustrated by Fig. [2j 

From the viewpoint of dynamic-order system, when investigating some multi-system prob- 
lems, we can establish the following generalized form of dynamic-order fractional dynamic 
systems 
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(3) 



The most important feature of the above generalized form is that the intricate interaction 
between dynamic systems has been represented by the differential order of each sub-system. 

In the following, we further illustrate the dynamic-order fractional dynamic system 
through two cases. 

Case 1. We firstly consider the relaxation system, which has been widely applied in energy 
dissipation, viscoelasticity and rheology, etc. 
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, |27[. To characterize the relaxation 



process, in which the relaxation pattern changes with control parameters or other variables, 
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the following variable-order fractional differential equation may be employed 

c Df Z) x(t) = -Bx(t) + f(t), < a(Z) < 1, 

(4) 

x(0) = 1, 

where B is the relaxation coefficient and Z is an independent variable. 

In engineering fields, fuzzy systems have extensive applications in control and system 
parameter calibration, and have become an efficient tool to exploit clear conclusion from 



uncertain information and incomplete data 



28N30|. Thereby, we consider the dynamic 



system fl4]) in which the variable-order is governed by the following fuzzy dynamic system. 
Here, we select the Takagi-Sugeno (T-S) fuzzy system as an example for illustration purpose. 



The T-S fuzzy system obeys the following rule 31] : 



Plant Rule v. If si(t) is nn and • • • and s p (t) is then 

y{t) = A i y (i),i = 1,2,... ,N, (5) 

where /i^- is the fuzzy set and N is the number of IF-THEN rules; si(t), . . . , s p (t) are the 
premise variables. Then the final output of T-S fuzzy system is inferred as follows: 

N 

y{t) = '52h i {y(t))Aiy(t), 

i=i (6) 
2/(0) = 1.0. 

In this illustrative case, we assume N = 2 and the fuzzy membership hi(y(t)) = | — 
\y{t),h 2 {y{t)) = \ + \y{t). 

If A\ = 0, Ai = — 1, then the exact solution of (J6j) is 



For numerical simulation, we assume that the relationship between the output of fuzzy 
system ([6]) and variable-order of fractional dynamic system (j2J) can be stated as follows 

a(Z) = a(t) = 1.0 -0.2y(t). (8) 

The evolution curve of variable-order a(t) and the corresponding numerical result of the 
relaxation system (j4j) are shown in Fig. [3l 
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FIG. 3: Left: The curve of the variable-order a(t) originated from the fuzzy system ©. Right: 
The relaxation curve of the system @ with the dynamic-order a(t) given in ([8]) and B = 1.0. 

From the observation of Fig. [31 the system flU) exhibits accelerating relaxation behavior 
with the variable-order (jSJ). Since fuzzy systems have been regarded as efficient tools to tackle 
real- world engineering problems, this case implies that, fuzzy systems may be an important 
method to characterize the evolution behavior of variable-order in fractional dynamic system. 

Case 2. We further consider the variable-order fractional diffusion system, in which 
the differential order is determined by a fractional order dynamic system. It means that 
the differential order in the considered fractional diffusion system relates to the output of 
another fractional dynamic system. It is a more general form to investigate the dynamic- 
order system. The considered variable-order fractional diffusion system is stated as 

< u(x, 0) = sin(x), x G [Q,L], (9) 
k u(0,t) = u(L,t) = 0,te [0,M], 

where K > is the generalized diffusion coefficient, u(x,t) is concentration, mass or other 
quantities of interest, q(x,t) is a source term and a(T) is a function of temperature. 
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FIG. 4: Left: The evolution curve of the temperature (jlOp with /3 = 0.9. In numerical simulation, 
the time step is 0.01. Right: The diffusion curve of the system Q with L = 1.0, M = 5.0, x = 
0.5, q(x,t) = and K = 0.1, the differential order is (fTTf) . 

We suppose the evolution process of temperature can be characterized by 

c D?T(t) = 0AT(t) + 10/(1.3* +1), 
T(0) = 10, 

where /3 G (0, 1] is the factional order derivative. Then, we assume the relationship between 
the differential order fl9]) and the temperature is 

a(T) = 0.8 + 0.005T. (11) 

It indicates that the diffusion process characterized by ((9]) is an accelerating subdiffusion 
process. The trajectory of a(t) and the diffusion curve of system (0) at x — 0.5 are drawn 
in Fig. HJ 

This example shows that the heat transfer process with result of temperature increasing, 
causes the accelerating behavior of considered diffusion process. In engineering situations, 
this model can offer an effective tool to explore physical mechanism of real-world diffusion 
related dynamic processes. 

Finally, we should note that, though we have made a first attempt to introduce the 
concept of dynamic-order fractional dynamic system, in which the multi-system interaction 
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is delivered by the differential order of each system. The further investigations on physical 
mechanism and application potentials are deserved. 
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